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ABSTRACT 

We obtained a second Chandra timing measurement of the 3.82 s pulsar CXOU J171405. 7— 381031 
in the supernova remnant (SNR) CTB 37B, which shows that it is spinning down rapidly. The 
average period derivative of (5.88 ± 0.08) x 10~^^ over the 1 year time span corresponds to a dipole 
magnetic field strength Bg — 4.8 x 10^'* G, well into the magnetar range. The spin-down power 
E = 4.2 X 10^'' erg s~^ is among the largest for magnetars, and the corresponding characteristic age 
Tc = P/2P — 1030 years is comparable to estimates of the age of the SNR. The period derivative 
enables us to recover probable pulsations in an ASCA observation taken in 1996, which yields a mean 
characteristic age of 860 years over the longer 13 year time span. The source is well detected up to 
10 keV, and its composite spectrum is typical of a magnetar. CTB 37B hosts HESS J1713— 381, the 
first TeV source that is coincident with a magnetar. While the TeV emission has been attributed 
to the SNR shell, it is possibly centrally peaked, and we hypothesize that this particularly young, 
energetic magnetar may contribute to the HESS source. We also searched for pulsations from another 
source in a HESS SNR, XMMU J173203.3-344518 in HESS J1731-347/G353.6-0.7, but could not 
confirm pulsations or long-term flux variability, making it more likely that this source is a weakly 
magnetized central compact object. 

Subject headings: ISM: individual (CTB 37B, G353.6-0.7) — pulsars: 

individual (CXOU J171405. 7-381031, XMMU J173203.3-344518) — stars: neu- 
tron 



1. INTRODUCTION 

Of the w 60 Galactic TeV sources"'^, the majority (43) 
are identified as either pulsar wind nebulae (PWNe) or 
supernova remnants (SNRs). The PWN class comprises 
« 27—30 of these, and it is likely that some of the unclas- 
sified sources will also prove to be PWNe. Most pulsars 
that are responsible for the PWNe detected at X-ray and 
TeV energies have spin-down power E = Att^IP/P^ > 
Iq3S gj.g g-i Figure[T]we display the distribution of E 
for rotation-powered pulsars powering TeV nebul ae. For 
reviews o f the se, sec Carri gan e t al. (2008), Galla nt et a,l.l 
(12001). iLemi'ere et all (I2008D. iHessels et al.l (120081) . 
iMattana et al.l (|2009( ). fc iKargaltsev fc Pavlov I |2010l ). 
In the standard model of nonthermal emission from 
PWNe, the X-rays are synchrotron from relativistic elec- 
trons and positrons, while the TeV 7-rays are inverse 
Compton scattered microwave background and other am- 
bient photons fr om the same population of relativis- 
tic pa rticles (e.g.. lde Jager fc Diannati-Ata'ill2008LlZhang| 
|2008| ). In contrast, there is scant evidence that any of 
the 17 known magnetars^, with periods in the range 
2—12 s, produce X-ray or TeV PWNe. Including both 
anomalous X-ray pulsars (AXPs) and soft gamma-ray 
repeaters (SGRs), the spin-down power of magnetars 
is < 10^^ erg s~^, as also plotted in Figure [H Be- 
cause of their larger magnetic field strengths and dif- 
ferent emission mechanisms from ordinary spin-powered 
pulsars, it is not predicted that magnetars accelerate par- 

VHE 7-ray Sky Map and Source Catalog, 



jhttp: //www. mpp mu.mpg.de/ rwagner/sources/index.html 

^ 5GR /AXP Online Catalog, 

|http : / / www . physics . mcgiU. ca/ ~ [ pulsar / magnetar /main, html 



tides to TeV ener gies (jBeloborodov fc Thomijo^ 120071 : 
i Thompsonl [2008l bl). except perhaps in their earliest 
sta ge of rapid spin-down (lAro ns 2003). 

iHalpern fc Gotthelj ()2010al hereafter Paper I) pre- 
sented a study of point X-ray sources in two SNRs 
detected at TeV energies by the HESS array of at- 
mospheric Cherenkov telescopes. The first of these, 
HESS J171 3-381 associated with CTB 37B, was dis- 
covered by lAharon ian et al. | (120061). Us ing a Chan- 
dra observation, lAh aronian et all (|2008aD then found 
the point X-ray source CXOU J171405. 7-381031 in 
CTB 37B, and considered that it could be a pulsar, 
albeit with an unusu ally soft, non-thermal spectrum. 
INakamua et al.l (j2009D analyzed Chandra and Suzaku 
spectra of CXOU J171405. 7-381031, noting possible 
variability in flux, which they took to be good evi- 
dence that it is an anomalous X-ray pulsar (AXP). 
In Paper I, we reported the discovery of pulsations 
from CXOU J171405. 7-381031 that verifles this conjec- 
ture. In Section 2 of this paper, we report a follow-up 
Chandra observation that measures the period deriva- 
tive of the pulsar, establishing its quantitative magne- 
tar properties.^ We raise in Section 3 the possibility 
that at an early stage, a magnetar such as the one in 
CTB 37B may produce a TeV PWN. In Section 4 we re- 
port on a new observation of the second compact X-ray 
source discussed in Paper 1, XMMU J173203.3-344518 
in HESS J1731-347/G353. 6-0.7, in which pulsations are 
not confirmed. 

While this paper was in preparation, the same result was an- 
nounced bv lSato et al.l l|2010l ) using an XMM-Newton observation. 
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TABLE 1 

Log of Timing Observations of CXOU J171405. 7-381031 in CTB 37B 



Instrument /Mode 


ObsID 


Date (UT) 


Date (MJD) 


Exp. (s) 


Counts (s ^) 


Period (s) 


ASCA GIS 

Chandra ACIS-S3/CC/F 
Chandra ACIS-S3/CC/F 
XMM-Newton EPIC pn 


54002030 
10113 
11233 
0606020101 


1996 Sep 12 

2009 Jan 25 

2010 Jan 30 
2010 Mar 17 


50,338.8 
54,856.3 
55,226.5 
55,272.2 


13,322 
30,146 
30,124 
40,264 


0.053 
0.106 
0.122 
0.264 


3.7954(1) 
3.823056(18) 
3.824936(17) 
3.825352(4)=^ 



From Sato ot aL (2010). 



2. X-RAY OBSERVATIONS OF 
CXOU J171405. 7-381031 IN 
HESS J1713-381/CTB 37B 

We obtained a new Chandra observation of 
CXOU J171405. 7-381031 on 2010 January 30 us- 
ing the same instrumental setup as our previous 
observation of 2009 January 25 (Paper I). Both data sets 
were acquired using the Advanced Camera for Imaging 
and Spectroscopy (ACIS) in continuous-clocking (CC) 
mode with the pulsar located on the ACIS-S3 CCD. This 
provides a time resolution of 2.85 ms and avoids any 
spectral pile-up. The photon arrival times in CC mode 
are adjusted in the standard processing to account for 
the known position of the pulsar, spacecraft dither, and 
SIM offset. Reduction and analysis used the standard 
software packages CIAO (v4.1) and CALDB (v3.5.0). 
A log of the Chandra timing observations is given in 
Table □ 

With a period derivative derived from the Chandra ob- 
servations (see below), we also searched for and likely 
detected pulsations in archival data tak en with the 
ASCA X-ray telescope (|Tanaka et al.lll994D . The pulsar 
was imaged with the Gas Imaging Spectrometer (CIS) 
on 1996 Se ptember 12 as part of the Galactic Survey 
(jSugizaki et al.H20dTI ). It fell on the very edge of the GIS 
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Fig. 1. — Spin-down power of 27 pulsars (hatched) whose 
PWNe are identified as TeV sources by HEGRA, HESS, VERI- 
TAS, MAGIC, and Milagro, in comparison with spin-down power 
of 17 magnetars (open squares) with measured period derivatives. 
The TeV PWNe are drawn from the reviews referenced in Section 
1, and magnetars from references in Footnote 2. Magnetars that 
have had SGR outbursts, including IE 1547.0—5408, are labeled in- 
dividually; the unlabeled squares are AXPs. The spin-down power 
of CXOU J171405. 7-381031 in CTB 37B is among the largest of 
magnetars, and falls in the range occupied exclusively by SGRs, 
although it is not known to be an SGR. 



field of view of sequence #54002030, 21' from the point- 
ing direction, partially in the area masked out by the 
standard processing. To maximize the sensitivity to the 
expected pulsar signal, we combined unmasked and un- 
filtered data taken in low- and high-bit-rate data modes, 
with timing resolution of 0.5 s and 0.0625 s, respectively. 
We obtained a total of 712 counts (0.5—10 keV) from the 
source in a 3' radius aperture during 13.3 ks of usable ex- 
posure from the 21 ks observation. 

2.1. Timing Analysis 

The photon arrival times were corrected to the so- 
lar system barycenter using the pulsar position obtained 
from the Chandra ACIS-I image. Timing analysis of the 
Chandra data was then performed in an identical manner 
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Fig. 2.— power spectra of CXOU J171405.7-381031 from 
the Chandra CC-mode observations separated by 1 year, and cor- 
responding folded light curves (inset). Normalized count rate is 
plotted. Phase is shifted arbitrarily. Note the highly significant 
decrease in frequency. 
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TABLE 2 

Chandra Measured Spin Parameters of 
CXOU J171405. 7-381031 in CTB 37B 



Parameter 



Value 



R.A. (J2000)'' 
Decl. (J2000)'' 
Epoch (MJD) 
Spin period, P 
Period derivative, P 
Range of dates (MJD) 
Surface dipole magnetic field. 
Spin-down power, E 
Characteristic age, Tc 



Bs 



17'U4™05=.74 

-38°10'30".9 

55226.5 

3.824936(17) s 

(5.88 ±0.08) X 10-" 

54856-55227 

4.8 X IQi* G 

4.2 X 10^* erg s"! 

1030 yr 



Measured from Chandra ACIS-I ObsID 6692. 
Chandra ACIS coordinate uncertainty is 0'.'6. 



Typical 



AP/At = 1.05(5)xl0" 



AP/At = 5.88(8)xl0" 



2009 



2010 



Year 



to Paper I. Photons from the central four source cohimns 
(2") in the range 1—5 keV, which contains most of the 
photons, were extracte d and search ed for pulsations us- 
ing the Rayleigh test (jStruttl [l88Cl[ ). also known as 
(|Buccheri et al.|[l983f) . The power-spectrum and folded 
pulse profiles from the two Chandra timing observations 
are shown in Figure [5] Although the eye may be drawn 
to subtle differences, a test on the difference profile 
indicates that the two light curves are not significantly 
different: = 1-44 for 20 degrees of freedom, corre- 
sponding to 9% chance probability. The pulsed fraction 
is fp w 0.38, defined as the fraction of counts above the 
minimum in the light curve, corrected for background. 
There is no evidence for energy dependence of the pulse 
profile. While it is possible that the true period is 7.6 s, 
with two peaks per rotation, there is no strong evidence 
to support this in the form of distinguishable peaks in 
the current data. 

The periods and Icr uncertainties for the two Chan- 
dra observations are listed in Table [TJ while the pe- 
riod derivative listed in Table [2] is calculated from the 
difference between the periods measured 1 year apart. 
The uncertainty on P is propagated from the statisti- 
cal errors in the individual period measurements. The 
P = (5.88 ± 0.08) X 10"" over the 1 year time span 
corresponds to a dipole magnetic field strength Bs — 
3.2 X 10i3(PP)i/2 = 4.8 X 10" G, well into the magnetar 
range. The spin-down power E = 4.2 x lO^'' erg s~^ is 
among the largest for magnetars, and the characteristic 
age Tc = P/2P = 1030 years is among the smallest. 

A third peri od measurement was obtained by 
iSato et al.l ()2010[ ) using XMM-Newton, only 47 days after 
our second Chandra timing observation. They reported 
P = 3.825352 ± 0.000004 s on 2010 March 17-18. As 
shown in Figure [31 this third point requires a significant 
increase in P to 1.05 x 10^^". Such variation is com- 
mon in magnetars, and in this case may indicate that 
the spin-down power was recently as large as « 7 x 10'^'* 
erg s~^. 

We recovered a fourth period measurement from the 
archival ASCA observation listed in Table [1] by extrap- 
olating the Chandra ephemeris. This detection shown 
in Figure SI while not significant in a blind search, has 
a false detection probability of only 1.1 x 10"^ when 
the search is restricted to a range of P up to twice the 
Chandra measured value. The implied period deriva- 
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Fig. 3. — (Top) Period measurements of 

CXOU J171405.7-381031 from the two Chandra CC-mode 
observations (filled circles) and the one XMM-Newton observation 
(Stato et al. 2010, open circle). Error bars are smaller than 
the size of the symbols. Period derivatives between consecutive 
measurements are indicated. A significant increase in P is seen. 
(Bottom) Average period derivative from 1996 to 2010 when 
including the probable ASCA CIS detection in 1996. The box 
encloses the range of the top panel. 



five relative to the 2009 Chandra measurement is P = 
7.08 X 10"^^, and the strength of the highly modulated 
signal is plausible. As shown in Figure |31 which in- 
cludes an extended range of period derivatives up to 
10 times the Chandra value, no other significant peak 
is found. Figure |3] shows the average P obtained in an 
unweighted fit to all four exsiting period measurements, 
P = 7.03 X 10^^° corresponding to Tc — 860 years, very 
close to the Chandra measured values, but perhaps a bet- 
ter representation of the long-term average. Keeping in 
mind that P is evidently variable, the ephemeris in Ta- 
ble [2] includes only the two Chandra points, which present 
the minimum observed value of P. 

2.2. Spectral Analysis 

Spectral data for the pulsar were extracted from the 
Chandra CC-mode observations using the five central 
source columns, corresponding to a diameter of 2."46 and 
« 95% of the point-source energy enclosed. The back- 
ground for the spectrum was obtained from the regions 
adjacent to the source. Background subtracted count 
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TABLE 3 

Chandra Spectral Fits for CXOU J171405. 7-381031 in 
CTB 37B 



Implied Period Derivative (10 



Fig. 4.— Period search of CXOU J171405. 7-381031 in the 
ASCA GIS observation of 1996 September 12. Period is converted 
to the P needed to extrapolate from the Chandra observation of 
2009. Inset: Folded light curve in the total (0.5-10 keV) energy 
band corresponding to the indicated highest peak in the power 
spectrum. Normalized count rate is plotted. 



rates are 0.106 ± 0.002 and 0.122 ± 0.002 s-\ respec- 
tively, for the two observations, indicating a likely change 
of « 15%. 

The standard spectral response matrices were gener- 
ated for the target location on the CCD using the CIAO 
tool psextract, with care taken to normalize the back- 
ground according to area. All spectra were corrected for 
the effects of charge-transfer inefficiency; however, the 
spectral gain is not calibrated well in CC mode. The 
resulting spectra were grouped with a minimum of 50 
counts p er bin and fitted using the XSPEC software 
IArnaudl (|l996l to power-law, two blackbody, and Comp- 
tonized blackbody models. The spectral fits are shown in 
Figure [SI and the best fitted parameters for these models 
are listed in Table [31 A single blackbody model does not 
yield a good fit, and is not shown. 

In paper I the ACIS-I/TE mode observation of 2007 in- 
dicated a rising spectrum above 6 keV, but the CC-mode 
spectra presented here do not clearly show that property. 
We suspect that the previous spectrum may have suffered 
systematic effects as the pulsar was dithered into the gap 
between CCDs. In addition, a small amount of spectral 
pileup in TE mode is avoided in CC mode. Nevertheless, 
the source is definitely detected up to 10 keV, and shows 
some evidence for an upturn there in the slightly brighter 
state of 2010. 

3. DISCUSSION 

3.1. CXOU .1171405.7-381031 as a Magnetar 

The Chandra CC-mode spectra of 
CXOU J171405. 7-381031 are fitted equally weU 
by power-law, two blackbody, and Comptonized black- 
body models. Of these, we disfavor the steep F « 4.2 
power-law model as being least plausible on physical 
grounds. In the thermal models, there is a blackbody 
component of temperature ~ 0.4 keV that covers ~ 10% 
of the surface area of the neutron star. A harder 
component is also needed that can be parameterized as 
either a second blackbody or a power-law tail, the latter 
approximating wh at is pro bably cyclotron resonance 
scattering (,Fcrnandez fc Thompson,2007.) using a simpli- 



Parameter 


2009 Jan 25 


2010 Jan 30 


Power-Law Model 


TVh (10^^ cm-^) 


5.46 ± 0.36 


5.36 ± 0.34 


F 


4.25 ± 0.23 


4.12 ±0.22 


F:r(2 - 10 keV)'' 


1.26 X 10-12 


1.48 X 10-12 


L^(2 - 10 keV)'' 


2.4 X 10^4 


2.7 X lO^-* 




1.06(56) 


1.07(67) 


Two Blackbody Model 

AT /i<i22 -2\ 01+0.44 ^ ni+U.tIV 



kTi (keV) 

Ri (km) 

kT2 (keV) 

i?2 (km) 

F^{2 - 10 keV)=- 

L^(2- 10 keV)'' 

r b 



■--0.27 
-+0.02 
-0.08 
^+3.1 
-•-0.9 
3+0.48 
-0.27 
3+0.32 
-0.13 



0.45: 

2.( 
1.13! 

o.ig!* ^ 

1.28 X 10-12 
1.8 X 10^'' 
4.2 X 10^4 
1.11(54) 



0.43^ 



Comptonized Blackbody Model 



-0.28 
n+0.03 
^-0.08 
o 9+3.9 
•-'•^-1.2 
1 -,0+0.38 
^•^^■^-0.23 

n r)(j + 0.28 

'J-^^-0.14 

1.50 X 10-12 
2.2 X 10^4 
5.4 X lO^-l 
1.07(65) 



TVh (1022 cm-2) 


4 12+"-''" 
^•1^-0.41 


4 1 l+U-«^ 

^■11-0. 34 


kT (keV) 


n oQ+0.08 


n qo + 0.08 


R (km) 


q q+3.4 
•^■'-'-1.6 


q c: + 9.0 


r 


o 04+0. 18 
■^•■^^-0.20 


■J-^^_0.15 


F^{2 - 10 keV)"^ 


1.32 X 10-12 


1.57 X 10-12 


L^(2- 10 keV)'' 


1.8 X lO^* 


2.2 X 10^" 




1.10(55) 


1.04(66) 



^ Absorbed flux in units of erg cm 2 s 1 . 

1" Unabsorbed luminosity in units of erg s^i, assuming d = 8 kpc. 



fied ana lytic expression in th e Comptonized blackbody 
model ([Halpern et al.l [20081 and references therein). 
The hot spot(s), and possibly anisotropic scattering in 
the magnetosphere, are responsible for the observed 
pulsations, although there is no obvious dependence of 
the pulse shape on energy. In either the two blackbody 
or the Comptonized blackbody model, the associated 
iVn « 4 X 10^^ cm^^ is compatible with its value in fits 
to the thermal emission from the SNR ([Nakamua et al.l 
[2009f ). which yield TVr = 3.5l^;^ x 10^^ ^^-2 rpj^jg 
provides some additional support for the association of 
CXOU J171405. 7-381031 with CTB 37B. 

Other properties of CXOU J171405. 7-381031 are con- 
sistent with those of AXPs: a factor of ~ 1.8 flux decrease 
is reported between the earlier S uzaku and Chandra ob- 
servations ([Nakamua et al.|[2009[ ). The spin-down power 
of CXOU J171405. 7-381031, S > 4.2 x lO^^ erg s'^, 
is among the highest for magnetars, exceeded only by 
SGR/AXP I E 1547.0-5408 (-B > 1 x 10^^ erg s'^; 
ICamilo et al] [2007. 2008), and equalled by SGR 162 7-41 
{E = 4.3 X 10^4 erg s'^; [Esposito et all [2001 ) and 
SGR 1806-20 when in its active state, - 5 x lO^^ 
erg s~^, although here we adopt a more typic al quiescent 
value of 8 X 10^3 erg s'^ ([Woods et al.[[2007[ ). The range 
of E in which CXOU J171405. 7-381031 falls (Figure [J) 
is occupied by SGRs rather than AXPs. By analogy with 
the history of IE 1547.0—5408 and other high-ii^ magne- 
tars that exhibit large fluctuations in P, we speculate 
that CXOU J171405. 7-381031 is prone to having SGR 
outbursts. It also shares with IE 1547.0—5408 the prop- 
erty that its spin-down power is greater than its X-ray 
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Fig. 5.— Spectra of CXOU J171405. 7-381031 in CTB 37B from the Chandra ACIS-S3 CC-mode observations of 2009 (left) and 2010 
(right )fitted to power- law (top), two blackbody (middle), and Comptonized blackbody models (bottom). Fitted parameters are listed in 
Tablets] 
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luminosity [L^ « 2 x 10^"' erg s~^), at least in quiescence. 
However, this does not change our assessment that the 
majority of its luminosity is generated by magnetic field 
decay and not by rotation, given its high surface dipole 
field strength, Bg — 4.8 x 10^'' G, X-ray spectral proper- 
ties, and variability that are typical of magnetars. 

3.2. Comparison of Age Estimates 

At radio wavelengths, CTB 3 7B is a shell « 10' 
in diameter (jKassim et al.l I1991D . Its distance, esti- 
mated from H I abs orption and t he Galactic rotation 
curve, is ~ 8 kpc (|Green I l2009f ). Estimates of the 
age of CTB 37B are quite varied, and it is of great 
interest to compare these with the spin-down age of 
CXOU J171405. 7-381031, even though characteristic 
ages of magnetars are not reliable because their period 
d erivatives are known to vary by a factor of a few. 

IClark fc StephensonI (|1975f) tentatively associated 
CTB 37B with the historical supernova of 393 AD, which 
would make it 1600 years old, almost twice the pul- 
sar's characteristic age of ~ 860 years. However, in 
the same region of the sky the more recently discovered 
SNR RXJ 1713.7-39 46 is a much bette r candidate fo r 
SN 393, as argued bvlWang et al.l (IT997h . IWaMi (l2006h . 
and ICassam-Chenai' et al.l ( 20041). At d ^ 1.3 kpc an d 
with TVh « 5 X 10^1 cm-2 (ICassam-Chenai et al.l[2Q0l . 
RXJ 1713.7—3946 is closer, and an order of magnitude 
less absorbed than CTB 37B, while the latter is too far 
away to have been visible for 8 months as the Chinese 
records indicate. Given this better alternative, we do not 
regard the birth of CXOU J171405. 7-381031 as having 
been established by historical records. Instead, we turn 
to estimates of the SNR age from analyses of its thermal 
X-ray spectrum. 

Using the Chandra data, lAharonian et al.l (|2008al) fit- 
ted a non-equilibrium ionization model and obtained a 
Sedov age of ~ 4900 yr for CTB 37B assuming electron- 
ion equilibrium, while noting that this condition is usu- 
ally not obtained in young SNRs. Both a higher ion 
temperature and efficient cosmic ray acceleration have 
the effect of decreasing the inferred age, pe r haps t o 
~ 2700 yr in the esti mate of lAharonian et aTl ()2008aD . 
iNakamua et al.l (l200tD then analyzed the Suzaku spec- 
trum of CTB 37B and estimated the age of the plasma 



as 650 



-2500 
-300 



years in the brightest part of the rem- 
nant, which also required a relatively low ISM density 
of riH ^ 0.4 cm~'^. In addition, they found a region of 
non-thermal power-law emission and a smaller ambient 
density in the southern part of the SNR. Thus, the SNR 
age estimate is broadly consistent with the pulsar's char- 
acteristic age, although the uncertainties in both values 
are large, and their applicability can be questioned. 

3.3. Magnetar/SNR Associations 

Three magnetars are securely associated with shell 
SNRs because of their central locations: IE 2259-1-586 
in CTB 109, IE 1841-045 in Kes 73, SCR 0526-66 in 
the LMC SNR N49 (Gaensler et al. 2001). To these we 
may tentatively add IE 1547. 0—5408 because of its poss i- 
ble radio shell G327.24-0.13 (IGelfand fc Gaensledl2007l) . 
and now CXOU J171405. 7-381031 in CTB 37B. In ad- 
dition to constraining ages for modelling of PWNe, such 
associations are valuable for inferring natal properties 



such as magnetar kick velocity, and possible contribu- 
tion of an assumed rapidly spinning proto-magnetar to 
the energetics of the SNR expansion (e.g., Vink & Kuiper 
2006). Theories for the birth of magnetars require short 
initial spin periods, of order a few milliseconds, in order 
to create their strong B-fields from a turbulent dynamo 
whose strengt h depends on the rotation ra te of the proto- 
neutron star (jThompson fc Duncanlll993l ). 

The radio emitting AXP XTE J1810-197 is the only 
magnetar with a directly measured proper motio n, cor- 
respo nding to Vt = 212 ± 35 (is. 5 km s~^ (Hclfan d et al.l 
I2007D . a value which is not unusual compare d to or- 
dinary young neutron stars. lAharonian et al.l (|2008aD 
noted that the displacement of CXOU J171405.7-381031 
from the center of the shell of CTB 37B would require 
a transverse velocity of ~ 1000 km s"^ for an age of 
5000 yr. The velocity would be unreasonably larger if 
the pulsar characteristic a ge of ^ 1000 years is ado pted 
as its true age. Howe ver, lAharonian et al.l ()2008aD and 
INakamua et al.l ()2009[ ) also noted that the brighter X-ray 
and radio emission on the eastern side of the remnant (see 
Figure 1 of Paper I) is probably due to a higher density 
there, implying that the pulsar may still be close to the 
explosion center, which is not the geometric center, re- 
ducing its inferred velocity. 

Because of this non-uniformity in shock radius and 
ISM density po, it is difficult to use the standard Sedov 
evolution 

2.026 Sot^ 



po 



as did iVink fc Kuipe^ ()2006D for other SNRs hosting 
magnetars, to estimate the energy Sq of the explosion. 
The radius, which enters to the fifth power, is uncer- 
tain by at least a fac tor of 2. Nevertheless , if we use 
the values adopted by INakamua et al.l ()2009( ) for the ra- 
dius, age, and ambient density in the brightest part of 
the remnant (their region 1), an unremarkable energy of 
£o ~ 1 X 10^° erg results. Energetic input to the SNR 
expansion from a presumed millisecond magnetar spin- 
nin g down is therefo re not evident, as was also concluded 
by ((Vink fc Kuiped |2006[) for other SNRs hosting mag- 
netars. 

3.4. Can the Magnetar Contribute to TeV Emission? 

HES S observations reported in lAharonian et al.l 
(2008a) locahzed the TeV emission to the center of the 
radio shell of CTB 37B, and found that the extent of 
the TeV source is compatible with either a centrally 
peaked Gaussian of a = 2.'6 ± 0.'8 or a shell of ra- 
dius 4'~6', consistent with the size of the radio remnant. 
The energy flux of HESS J1713-381 from 0.2-10 TeV is 
« 4.2 X 10~^^ erg cm^^ s^^, corresponding to luminosity 
Lj = 3.2 X lO'^* dg erg s~^, which we note is almost equal 
to the spin-down power of CXOU J171405. 7-381031, 
E = 4.2 X lO'^'* erg s^^. Previous authors have argued 
that particles accelerated at the SNR shock are respon- 
sible for the TeV emission fr om HESS J1713-381 via 
pion decay (jAharonian et al.] [200 8a) or multi-zone in- 
verse Compton scatter ing of the microwave background 
(|Nakamua et al.l 120091 ). As it is not clear if the TeV 
source has a pure shell-like morphology, it cannot be 
ruled out that the pulsar makes some contribution to the 
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TeV flux, possible via inverse Compton scattering from 
a PWN that is confined within the SNR shell. 

Since the TeV luminosity of HESS J1713-381 is 
almost equal to the present spin-down power of 
CXOU J171405. 7-381031, a scenario for a pulsar- 
powered TeV nebula would necessarily involve particles 
injected at an earlier time when its spin-down power was 
larger. However, unlike ordinary pulsars, it is not cer- 
tain that magnetars can accelerate particles to TeV en- 
ergies. Magnetar models that involve strong currents on 
closed, twiste d magnetic field lines develop vo ltages of 
only ~ 10^ V (|Beloborodov fc ThompswH [20071 ). On the 
other hand, it is not excluded that the ordinary pulsar 
mechanism operates on open magnetic field lines of mag- 
netars, which may generate particle-dominated winds 
that become shocked PWNe, accelerating a fraction of 
the particles to TeV energies. Evidence for this possibil- 
ity comes from the three magnetars that are also tran- 
sient radio pulsars (jCamilo et al.ll2006ll2007l : ILevin et al.l 
[2OT0h . Even so, we c aution that models for radio emis- 
sion from magnetars (jThompson 2008a, b) produce aver- 
age particle energies of only 3 GeV even on open field 
lines, because 7-rays make pairs more easily in the strong 
magnetic field. 

Nevertheless, if high-energy electrons are generated 
in a shocked PWN in an early stage and escape the 
high B-field region, they can emit via inverse Comp- 
ton scattering into the TeV band for an extended time 
(|de Jager fc Diannati-Ata^ [2008). Electron lifetimes 
against synchrotron and inverse Compton losses are 

4200 yr / B 

and 

_ 8 X lO'' yr 

respectively, for an electron producing a photon of en- 
ergy ^^..^(TeV) by inverse Compton scattering of the mi- 
crowave background. The limiting lifetime of such a relic 
nebula is then determined by tg as long as i? > 2 fiG. 

The spin-down power of CXOU J171405. 7-381031 was 
certainly large enough to generate a TeV PWN during 
the history of CTB 37B, as long as it was born with a 
period significantly shorter than its present value, since 
E = 9x 1036 p-4 Ibs/5 X lO^"^ G)2 erg s'^. The current 
TeV luminosity of HESS J1713-381, « 3.2x10^4 ergs'^ 
if radiated for 1000 yr, amounts to ~ 1 x 10*^ erg, which 
is < 10"'^ of the initial rotational energy of the pulsar 
for Pq < 0.14 s. Thus, we consider an inverse Comp- 
ton PWN to be a plausible channel for TeV emission 
from CXOU J171405. 7-381031 because it requires only 
a small fraction of the rotational energy to be deposited 
in very high-energy electrons. 

A possible prototype of such systems is 
PSR J1846-0258, th e "transitional" - 700 yr old pulsar 
(|Gotthelf et al.|[20i00 l) associated with HESS J1846-029 
in the shell-like SNR Kes 75 (Djannati-Atai 2008). This 
0.32 s pulsar has a dipole field that approaches magnetar 
strength {B^ = 4.9 x 10" G) and displays AXP-like 
bursts (Gavriil et al. 2008); it is therefore likely of an 
intermediate class connecting the rotation-powered and 
the magnetar pulsars. It is also one of the most energetic 



pulsars; its E — 8.1 x 10^^ erg s ^ is easily sufficient to 
power its X-ray PWN and the compact T eV emission 
of HESS J1846-029 (|Diannati-Atai [20081) . However, 
Kes 75 differs from CTB 37B in that it contains a 
powerful X-ray PWN, while there is no evidence of 
an X-ray PWN in CTB 37B. Although these SNRs 
are of comparable age, any X-ray PWN in CTB 37B 
may have evolved a factor of 10 faster than the one in 
Kes 75, based on the scaling of pulsar E with true age 
T such that T oc B~^ for a given E. The lifetime of 
a hypothetical rotation-powered X-ray PWN around 
a magnetar m u st be relatively brief. Interestingly, 
IVink fc Bambal (|2010D claim to have detected the first 
X-ray PWN around a magnetar in an archival Chandra 
observation of IE 1547.0—5408, the magnetar of highest 
E. Although this would be an important precedent, 
we hesitate to interpret it as such because we have not 
clearly confirmed its existence in those data. 

Evidence to associate any other magnetar with TeV 
emission is sparse. The AXP IE 1841-045 in SNR 
Kes 73 is located at t he edge of HESS J1841-055 
(jAharonian et al.l I2008bf) . but this is not a likely asso- 
ciation because the TeV emission is w 1° in diameter 
and may co nceivably hav e more than one source, as con- 
sidered by iSguera et al.l (|2009D . Moreover, it can be 
expected that any high-energy particles are still con- 
tained within the SNR shell, which is much smaller 
than the TeV source. Recently, an extended TeV source 
that surrounds the massive, young star cl uster West- 
erlund 1 has also been reported by HESS ()Ohm et al.l 
|2009( ). It is possible that this source is associated 
with a young pulsar born in th e cluster, the traii sient 
AXP CXO J164710. 2-455216 (IMuno et al.l 120061) . or 
its supernova remnant. Notably absent among the re- 
ported HESS sources is IE 1547.0-5408, the magnetar 
of highest E and a small characteristic age similar to 
CXOU J171405. 7-381031. It would be interesting to 
target IE 1547.0-5408 for a deep observation by HESS 
as a second test of whether a young magnetar can pro- 
duce a TeV source. 

4. X-RAY TIMING OBSERVATION OF 
XMMU J173203.3-344518IN HESS J1731-347 

In Paper I, we also presented an analysis of the 
point source XMMU J173203.3-344518 in the SNR 
G353.6— 0.7 that ha s been identified with the TeV source 
HESS J1731-347 (lAharonian et al.l I2008bl: ITian et al.l 
[2OO&) . We found only marginal evidence for a 
1 s period in an XMM-Newton observation, which, 
in combination with the absence of strong evidence 
for variability did not allow us to definitively clas- 
sify XMMU J173203. 3-344518. Both a magne- 
tar and a weakly magnetized central compact ob- 
ject (CCO) remained possibilities. The spectrum of 
XMMU J173203. 3-344518 was best fitted by a two 
blackbody model, which is common to members of both 
classes. 

We obtained a new Chandra observation of 
XMMU J173203.3-344518 on 2010 May 18 using 
ACIS-S in CC-mode for 40 ks. Using the same analysis 
techniques as described in Section 2 for CC-mode data, 
our period search of XMMU J173203.3-344518 did 
not reveal the candidate signal at 1 s, nor any other 
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TABLE 4 
Chandra Observation of 
XMMU J173203. 3-344518 in G353.6-0.7 



XMMU J173203.3-34451t 
2010 MAY 18 



Parameter 



2010 May 18 



ObsID 
Mode 

Exp. time (s) 
Counts (s~^) 
R.A. (J2000)'' 
Peel. (J2000)'' 



11234 
ACIS-S3/CC/F 
39,920 
0.334 
17''32™03=.40 
-34°45'16".77 



A^H (10^^ cm-2) 
fcTi (keV) 
Ri (km) 
kT2 (keV) 
R2 (km) 

F^(0.5 - 10 keV)'' 
L^{0.^ - 10 keV)'= 

^bol 



Two Blackbody Model 



2.13^ 



-0.14 

0.36+°04 
2.1^ 



-0.05 
hl.3 
-0.4 

0.61 ±0.07 

^+0.32 
-0.29 



0.44^ 



2.74 X 10- 



-12 



1.26 X 10^* 
1.34 X lO^"* 
1.066(102) 



Measured from Chandra ACIS-I ObsID 
9139. Typical Chandra ACIS coordinate un- 
certainty is 0'.'6. 

^ Absorbed flux in units of erg cm~^ . 

Unabsorbed luminosity in units of erg s~^, 
assuming d = 3.2 kpc. 



significant period. The 95% confidence upper limit on 
the pulsed fraction of a sinusoidal signal is 8.6% for any 
period down to 10 ms. If restricted to periods > 1 s, the 
corresponding limit is 7.6%. We also performed searches 
in restricted energy bands, in case there is a pulsation 
with an energy-dependent p hase reversal similar t o PSR 
J0821-4300 in Puppis A (IGotthelf fc Halpernl l2009l) . 
but no such signal was found. 

We fitted the spectrum to a two blackbody model, 
which was the best fitting model to previous observa- 
tions of XMMU J173203. 3-344518 in Paper I. The re- 
sult in shown in Figure [6l the data is grouped with 100 
counts per bin. Table 2] gives the parameters of the two 
blackbody fit. These are consistent with the previous ob- 
servations; in particular, the 0.5—10 keV flux is identical 
to that of an XMM-Newton observation from the year 
2007 (Paper I) . The absence of variability and pulsations 
tends to favor a CCO; seve ral of the latter have very 
low limits on pulsed fraction ([Halpern k, Gotth"el3l2010bl 
and references therein). If so, XMMU J173203.3-344518 
would be the most luminous CCO, with ~ \ y. lO'^^ 
erg s~^. Recently, HESS was able to resolve the TeV 
emission of HESS J 173 1—347 into a shell t hat matches 
the size of the radio remnant G353.6— 0.7 (|Acero et al.l 
|2P09). This would tend to rule out the central source 
XMMU J173203.3-344518 as a source of the TeV emis- 
sion. 

5. CONCLUSIONS 

In Paper I, we reported the discovery of the 3.82 s 
pulsar CXOU J171405. 7-381031 in CTB 37B, which 
is coincident with the TeV source HESS J1713-381. 
Additional X-ray timing measurements reported 
here for the first time reveal the spin-down rate. 
CXOU J171405.7-381031 is one of the most energetic 
magnetars in terms of spin-down power, and has a 
characteristic timing age of ~ 1000 years or less. Its 




TWO BLACKBODY MODEL 



4^ 



ENERGY (KEV) 

Fig. 6.— Spectrum of XMMU J173203.3-344518 in G353.6-0.7 
from the Chandra ACIS-S3 CC-mode observation of 2010 May 18 
fitted with two blackbodies. Fitted parameters are listed in Table|4] 



timing age is roughly consistent with other estimates of 
the age of the SNR. In addition, there is some evidence 
for long-term variability by a factor of ^ 2 in spin-down 
rate, which recently is found at its highest measured 
value. Its E falls in the range occupied by SGRs, which 
may portend a future SGR outburst. The off-center 
location of the pulsar inside the SNR shell is probably a 
result of the inhomogeneous medium in which the shell 
expanded, rather than indicating an unprecedented high 
natal kick velocity. 

There are considerable theoretical obstacles to getting 
particles of TeV energy from a magnetar, as well as hav- 
ing them emit primarily TeV 7-rays. As well, there is lit- 
tle observational evidence for such particles in the form 
of an X-ray synchrotron nebula. Nevertheless, the excep- 
tional youth and high E of CXOU J171405. 7-381031 
lead us to entertain the possibility that it generated a 
relic PWN that is now emitting primarily via inverse 
Compton scattering at TeV energies. Future imaging at- 
mospheric Cherenkov observations to try to resolve the 
shell of CTB 37B from the pulsar could test this hypoth- 
esis. 

Turning to another compact source in a 
HESS SNR, we obtained a second X-ray tim- 
ing observation of XMMU J173203.3-344518 in 
HESS J1731-347/G353.6-0.7, but did not succeed in 
discovering a period. Its nature remains ambiguous 
because its composite spectrum, best fitted with two 
blackbodies, is characteristic of either AXPs or CCOs. 
In comparison with three previous X-ray observa- 
tions of this source, the absence of obvious variability 
so far is consistent with XMMU J173203.3-344518 
being a CCO (anti-magnetar). Deeper observations 
will be necessary to discover its timing properties, 
which would determine conclusively whether it is a 
magnetar or an anti-magnetar. The present upper 
limit on pulsed fraction is comparable to the weakest 
pulsations observ ed in magnetars inc luding, at times, 
IE 1547.0-5408 ()Halpern et a l.l [200l . and similar to 



the up per limits on several CCOs (|Halpern fc Gotthell 
l2010bl) for which pulsations have not been detected. 
Even if pulsations are absent, continuing monitoring 
of XMMU J173203.3-344518 for variability may also 
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resolve its nature, since only a magnetar would be 
variable. 
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